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Abstract 

We impose the periodicity conditions corresponding to the Matsubara 
formahsm for Thermal Field Theory as constraints in the imaginary 
time path integral. These constraints are introduced by means of time- 
independent auxiliary fields which, by integration of the original vari- 
ables, become dynamical fields in the resulting 'dual' representation for 
the theory. This alternative representation has the appealing property 
of involving fields which live in one dimension less than the original 
ones, with a quantum partition function whose integration measure is 
identical to the one of its classical counterpart, albeit with a different 
(spatially nonlocal) action. 

1 Introduction 

Quantum Field Theory (QFT) models with constrained configuration spaces 
naturally arise within the context of modern applications, particularly gauge 
invariant systems [T|. 

More recently, a formulation involving constraints has been also applied 
to deal with different kinds of problems, namely, the static and dynamical 
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Casimir effects [21 [3l 0], and to derive tfie overlap Dirac operator [5], [6] in a 
simpler way [?]. There, fields used to impose the constraints lead, after inte- 
grating out the original variables, to an effective model where the dynamical 
fields live on the constrained surface. Since the codimension of the bound- 
aries is usually equal to 1, the effective model is defined in one dimension 
less than the original one. 

Here, we present the extension of that kind of approach to a less natural 
realm, that of QFT at finite temperature (T > 0), to deal with the periodicity 
constraints in the imaginary time. The foundations of QFT at T > were 
laid down quite a long time ago [HI [9]. The original approach to this topic, 
the now called Matsubara (or 'imaginary-time') formalism has been very 
successful indeed in allowing for the evaluation of thermal effects in QFT, 
both in High Energy [TUI and Condensed Matter Physics. 

It allowed, for example, to study the new phenomena that emerge when 
using a Statistical Mechanics description for quantum relativistic systems. It 
also provided a convenient way to naturally extend the notion of Abelian and 
non Abelian gauge fields, studying its consequences for particle physics [TOl 
[m [H], in the T > context. 

A fundamental property introduced by this formalism is the imaginary- 
time periodicity (antiperiodicity) conditions for the bosonic (fermionic) field 
configurations in the path integral. That may be clearly seen already at the 
level of the partition function, 2{j3), for a system at a temperature T = 
with a Hamiltonian H: 

Z{P) = Tr(e-'^^) . (1) 

Assuming first, for the sake of simplicity, that there is only one (bosonic) 
degree of freedom, described by a coordinate q, the expression above may be 
written more explicitly as follows: 

Z{(3) = j dq {qle-f'^lq) = J dq {q, -z(3\q,0) , (2) 

where \q,t) denotes the usual 'rotating basis' elements, which here appear 
evaluated at imaginary values of t. Then the standard path integral construc- 
tion for the transition amplitude between different times may be applied, to 
obtain the partition function in the Matsubara formalism: 

Z{P) = I pppge^'''"['^'^-^M , (3) 



'q(0)=q(/3) 

es phase 

that g(0) = g(/3), while the pij) paths have free boundary conditions 



where the measure includes phase-space paths qij), pij) (r e [O,/?]) such 

0. 



^Note, however, that a sHghtly more symmetric form for those conditions in the path 
integral for Z{(3) could be used [l]. 
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When applied to a bosonic field theory in d + 1 spacetime dimensions, 
this procedure leads to field paths which are periodic in the imaginary time, 
while the canonical momentum ones are, again, unrestricted. Moreover, when 
the Hamiltonian is quadratic in the canonical momentum, integration of this 
variable yields a model where the dynamical field is defined on S^xR^, where 
the radius of is proportional to the inverse temperature, (3. In Fourier 
space, the corresponding frequencies become the usual discrete Matsubara 
frequencies. 

A characteristic feature of the Matsubara formalism (shared with the 
real-time formulation) is that the introduction of a time dependence for the 
fields seems to be unavoidable, even if one limits oneself to the calculation of 
time independent objects. 

With the aim of constructing a new representation where only static fields 
are involved, we shall introduce here an alternative way of dealing with T > 
QFT calculations. The procedure is inspired by a recent paper in which a 
constrained functional integral approach is used to implement the eff'ect of 
fluctuating boundaries in the Casimir effect [2]. In the present context, this 
allows one to introduce the periodicity conditions by means of Lagrange 
multipliers (d-dimensional when the field lives in d + 1 dimensions). Then 
the original fields can be integrated, what leaves a functional depending only 
on the d-dimensional Lagrange multipliers. 

This paper is organized as follows: in section [2] we introduce the method, 
using the harmonic oscillator as a convenient framework. In section [3] we deal 
with the real scalar field, and in [4] a Dirac field is considered. In section [5] we 
present our conclusions. 



2 The method 

2.1 The periodicity constraint 

Let us see, again within the context of a system with a single degree of 
freedom, how the thermal partition function may be obtained by imposing 
appropriate constraints to the path integral for Zq, the (zero temperature) 
vacuum persistence amplitude. For reasons that will become clear below, we 
start from its phase-space path integral: 

Zo = y PpPge-'^['?(")'f(")l , (4) 

where S is the first-order action, S = dr with C = —ipq + H{p,q), 
and H denotes the Hamiltonian, assumed to be of the form: H{p, q) = 
T{p) + V{q). 
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Of course, Zq is the limit of an imaginary-time transition amplitude, 



Zo = lim {qo,-iT\qo,iT) 

T — >-\-oo 

= lim y |(go|n)pe-2^^" = lim Kgo|0)pe-2^^« 

n 

where we have introduced |n), the eigenstates of H, H\n) = En\n), and go, 
the asymptotic value for qo at T ^ ±00 (usually, go = 0). Eq is the energy 
of |0), the ground state. 

Let us now see how one can write an alternative expression for Z{P), by 
starting from the vacuum transition amplitude, Zq, and imposing the appro- 
priate constraints on the paths. To that end, we first use the superposition 
principle, introducing decompositions of the identity at the imaginary times 
corresponding to r = and r = P, so that we may write Zq in the equivalent 
way: 

Zo = lim / rfg2(igi (go,-iT|g2,-2/9) (g2,-2/3|gi,0) (gi,0|go,iT) , (6) 

T->oo J 

or, in a path integral representation, 

rqiT)=qo 



Zq = lim dq2dqi / VpVqe'Jp'^^^ 

J JqW=q2 

rq{P)=q2 „ rq{o)=qi „ 

X / VpVqe-^o^^^ / VpVqe~^-T^^^ . (7) 

Jq{0)=qi Jq{--T)=qo 

The representation above is quite useful in order to understand which is the 
correct way to impose the constraints, to obtain Z{P). In short, to reproduce 
Z{P) we have to impose periodicity constraints for both phase space variables. 
Indeed, let us introduce an object Zs{P) that results from imposing those 
constraints on the Zq path integral, and extracting a Zq factor: 

^ ... _ / VpVq 6{qi(3) - g(0)) - p(0)) 

= JVpVqe-^ • 

Then, the use of the superposition principle yields: 

J VpVq S{q{p) - g(0))5(p(/3) - p{0)) e^^ = ^lim j dp^dq^ 

(go, -iT\pi, -i/3) {pi, -z/5|gi, -i(3) (gi, -i/3\qi, 0) (gi, 0|pi, 0) (pi, 0|go, iT) 

9) 
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VpVq6{q{(3) -q{0))6{p{(3) -p{0))e~^ = lim e-^^^^r-/?) 



< j ^^(go|0)(Obi) (gi,-2/9|gi,0)(pi|0)(0|go) 
= lim e-^°(2T-/3)|^^^|Q^|2 fdq,{q,,-if3\q,,0). 

T^oo J 

= ZoX e^^°Z(/5) = Zox Tr [e'^^^"^")] . 




(10) 



Then we conclude that 




(11) 



: H := H - Ei 



(12) 



The conclusion is that, by imposing periodicity on both phase space vari- 
ables, and discarding /5-independent factors (since they would be canceled 
by the normalization constant) we obtain Zs{(3), the partition function cor- 
responding to the original Hamiltonian, the ground state energy redefined to 
zero. The subtraction of the vacuum energy is usually irrelevant (except in 
some exceptional situations), as it is wiped out when taking derivatives of 
the free energy to calculate physical quantities. 

Note that the introduction of periodicity constraints for both variables 
is not in contradiction with the usual representation, ([3l), where they only 
apply to q, since they corresponds to different sets of paths. Indeed, in our 
approach the new constraints are crucial in order to get rid of the unwelcome 
factors coming from paths which are outside of the [0, j3] interval (which are 
absent from the standard path integral). 

We conclude this derivation of the boundary conditions by showing ex- 
plicitly why the usual procedure of introducing a periodicity constraint for 
just the coordinate g(r) would not be sufficient. Indeed, we can see that 




)]} 




^""'yif^-f'^''^ - /<^*i(oi,.)n,..-^^i*.o> . (14) 
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Figure 1: Representation of the compactification mechanism 

where we cannot extract a 2{P) factor, due to the presence of the squared 
vacuum wave function inside the integral. It is not difficult to realize that 
that factor, whose entanglement makes it impossible to extract the partition 
function, is due to contributions from paths outside of the (0, P] interval. 

Summarizing, we have shown that the proper way to extract the parti- 
tion function from the T = partition function Zq, is to impose periodicity 
constraints for both the coordinate and its canonical momentum, a proce- 
dure that yields a Zq factor times the thermal partition function, Zg{(3). 
We present, in Figure 1, a pictorial representation of this 'compactification' 
mechanism. 

2.2 Auxiliary fields 

Let us now see how the use of auxiliary fields to exponentiate the constraints 
leads naturally to an alternative representation. The two (5-functions require 
the introduction of two auxiliary fields, and ^2- which are just real (time- 
independent) variables in this case. Using the notation Qi = q and Q2 = p, 
we have 



Using this representation for the constraints and interchanging the order of 
integration for the multipliers and the phase space variables, the resulting 




(15) 



6 



expression for Zs{i3) may be written as follows: 

where = Zq, and we have introduced the notation: 

Ja{r)^UHr-f3)-6{r)] . (17) 
The phase-space measure has been written in terms of Q: 

VQ^ n (18) 

— oo<r<oo 

For the particular case of a harmonic oscillator with unit mass and frequency 
cj, we have 
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hoo 



SiQ) = SoiQ) = - / drQair))CabQaiT) , (19) 
where ICab are the elements of the 2x2 operator matrix IC, given by: 

£^(4f). 

Thus the integral over Q is a Gaussian; it may therefore be written as follows: 

Zs{l3) = (det£)~^ j ^ e-s?-^^^-?" , (21) 

with 

M = fi(o+) +fi(o_) - n{p) - n{-p) , (22) 

where f2(r) denotes the inverse of the operator IC of (l20|l : namely, 

jCac^cbir) = SabSir) (23) 

where the fiat's denote the matrix elements of Q. 

The explicit form of this object may be easily found to be the following: 

' fsgn(r) \ ^^|,| 



(sgn = sign function) 
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Equation (|2^ can be used in (j22l) . to see that: 



M = [n{0+) +fi(0_)] (1 -e"^^ 

cj-i 




n^H + 1)"' , (25) 



where 

nsico) = (e^- - l)-i (26) 

is the Bose-Einstein distribution function (with the zero of energy set at the 
ground state). 

Finally, note that Af exactly cancels the ( det /C) ^ factor, and thus we 
arrive to a sort of 'dual' description for the partition function, as an integral 
over the variables: 

Z,{(3) = / — e" ^msM+u . (27) 
J 27r 

This integral is over two real variables ^a, which are 0-dimensional fields, 
one dimension less than the + 1 dimensional original theory. To interpret 
this integral we may compare it with the one corresponding to the classical 
statistical mechanics version of this system. To that end, we evaluate the 
partition function in the classical (high-temperature) limit. In that limit, we 
approximate the integrand accordingly to see that Zs{P) becomes: 



(28) 



where: 

H{^u^2) ^ l{ei + u;'e2)- (29) 

We see that (l28l) corresponds exactly to the classical partition function for a 
harmonic oscillator, when the identifications: C,i = p (classical momentum), 
and ^2 = 1 (classical coordinate) are made 

ZsW) - /^e-^^(-^ + -V) (;5«l). (30) 

On the other hand, had the exact form of the integral been kept (no ap- 
proximation), we could still have written an expression similar to the classical 
partition function, albeit with an 'effective Hamiltonian' iJe//(^i, ^2): 

MP) = I ^e-^^^/Ha,6)^ (31) 
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where: 

^e//(ei,6) = ^ {nsiuj) + 1)-' {u-'ei + ^ ^2') • (32) 

This shows that the quantum partition function may also be written as a 
classical one, by using a /3-dependent Hamiltonian, which of course tends to 
its classical counterpart in the high-temperature limit. 

By integrating out the auxiliary fields in the (exact) expression for the 
partition function (|27|) . we obtain: 

Z,{(3) = riBiu) + 1 = ■ (33) 

which is the correct result. In what follows, to simplify the notation, we 
shall omit writing the 's' subscript in Z{j3), assuming implicitly that one is 
dealing with the normal-ordered Hamiltonian. 

An important fact that has emerged from an analysis of the classical 
(high-temperature) limit: the auxiliary fields do have a physical interpreta- 
tion. The multiplier associated to the periodicity condition for q plays the 
role of a classical momentum, while the one corresponding to the periodic- 
ity for the momentum becomes a generalization of the classical coordinate. 
The same interpretation might also be retained far from the classical limit, 
but then the Hamiltonian departs from the classical one, receiving quantum 
corrections. 

This representation is also valid for interacting theories. To that effect, 
note that, even when the action S is not quadratic, we may still give a 
formal expression for the alternative representation. Indeed, denoting by 
Z{J) the zero-temperature generating functional of correlation functions of 
the canonical variables: 

Z{J) = j VQ e~^^^^^^-'^'^^'^''^^'>^^^^^ (34) 

and by W( J) the corresponding functional for connected ones, we see that 

= [Z{Or'l ^ exp{W[^j(r)]} , (35) 

where, with our normalization conventions, Z{0) = Zq (the vacuum func- 
tional for the interacting case). 

Thus, a possible way to derive the effective Hamiltonian in the interacting 
case is to obtain first W[J], and then to replace the (arbitrary) source J(r) 
by «j(r), where j(r) is the function of the auxiliary field defined in (fT71) . Of 
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course, W cannot be obtained exactly, except in very special cases. Other- 
wise, a suitable perturbative expansion can be used. In any case, W can be 
functionally expanded in powers of the source J(t): 

oo „ 

W[J] = W[0] + ^- / wi:)(ri,...,r„)Ja,(ri)...J.„(r„) (36) 

„=1 -^n,.-,T„ 

where each coefficient W'-"'' is the n-point connected correlation function. 
The expansion above immediately yields an expansion for -f^e// in powers 
of the auxiliary fields. It is important to stress that this expansion is not 
necessarily a perturbative expansion. Indeed, the strength of each term is 
controlled by W^"-*, which could even be exact (non-perturbative) in a cou- 
pling constant. To fix ideas, let us see what happens when one keeps only 
up to the n = 4 term, assuming also that there is Qa —Qa symmetry in 
S. Then, we first see that the W[0] is cancelled by the J\f factor, and on the 
other hand we obtain 



(27r) 



= I 7^e-^^^//(«-«^), (37) 



where 



~ 4^5 / ^Si2a3aA^l,^2,r3,T4)jaAn)ja2ir2)jasiT3)ja^{n) 

+ ... (38) 

Using the explicit form of jai'T) in terms of the auxiliary fields, we see that: 



Heff = Hf^f + H'fff + • • • (39) 



where 



^eff - 2^ab^a^b 
1 



2 



rr{2k) 
^eff 



where the explicit forms of the coefficients A^^^'^) in terms of W^^'^-* may be 
found, after some algebra. For example Ai^"^^ is a diagonal matrix: 

- ( I ) («) 



10 



where 

ca = \ / - (1 - e-^) Waail^) (42) 

P J 

(where the tilde denotes Fourier transform). It is immediate to realize that 
Cl plays the role of an effective coefficient for the kinetic term (oc p^) in the 
effective Hamiltonian, while C2 does introduce an effective quadratic poten- 
tial. Note that they will, in general, depend on P, u, and on any additional 
coupling constant the system may have. For the harmonic oscillator case we 
have the rather simple form: 

1 

Cl - 



uj{nB{uj) + 1) 

C2 = . (43) 

nB{oj) + 1 

The quartic term involves AA^^\ which may be written in terms of the con- 
nected 4-point function: 

1 

- 6 >vi;L(/?, /5, 0, 0) + 4 >v1:L(/?, 0, 0, 0) - w1:L(0, 0, 0, 0)1 (44) 

J sym 

where the sym suffix denotes symmetrization under simultaneous interchange 
of time arguments and discrete indices. Of course, this expression could also 
be written in Fourier space; we shall leave the analog of this term for the 
case of the real scalar field. 



->vSL(0,0,0,0) + 4>v1:L(A/3,A0) 



2.3 Generating functional 

Now we proceed to the calculation of thermal correlation functions within 
the approach that we developed for the calculation of the partition function. 
To that end, we shall introduce the generating functional of correlation func- 
tions, to be denoted by Zs{[3, J), and Ws(/9, J) = \nZs{P, J) the generating 
functional of connected correlation functions. With these conventions, the 
connected correlation functions are given by: 

From the previous subsection, we know that path integral expression for 
the generating functional shall be: 
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where S{Q) denotes the first-order form of the action. 

Let us evaluate ^s(/3, J) explicitly for the case of the harmonic oscillator, 
where iS is a quadratic form. The result of performing the Gaussian integral 
over Q, may be written in this follows: 

Using the explicit form for ja, we see that the expression above is equiv- 
alent to: ^ 

Z,(/3,J) = Z,{J) [ lie"^«»^'"'^' + *^'^^'^ (48) 
J 27r 

where Zs{J) is the zero-temperature generating functional, 

Zs{J) = '^'"^ '^'^'^ MTi)^ab{Ti,T2)M'r2) (^49^ 

and 

/+CO 
dT [(/C-^),,(/3, r) - (/C-i)„,(0, r)] Mr) . (50) 



Integrating out the auxiliary fields, and recalling the results of the previous 
subsection, 

ZsiP, J) = Z{(3) Zs{J) e"^^" [^i-'UN, (51) 

where the Na are the functionals of Ja defined in (1501 ). 

Neglecting a source-independent term (irrelevant for the calculation of 
correlation functions), the W generating functional will have the structure, 

>V.(/3,J) = ^JdnJ dT2JaiTi)Gab{Tl,T2)MT2) (52) 

where Gab denotes the thermal correlation function: 

Gabin, T2) = Gabin - T-2) = {Qain)QbiT2)) , (53) 
whose structure we shall now write more explicitly. We first note that: 

Gabin, T2) = Gi;)(n, r^) - t/if (n, r^) (54) 

where G^^^ is the zero-temperature correlation function while U^^^ denotes a 
temperature-dependent piece. 

We can write a more explicit form for the two terms that enter into the 
expression above for the correlation function. Indeed, for G^f^ we have: 

^ ^ ( _L ) (55) 
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while for f/^^^ we may use the following matrix representation: 

[/(/^) = [/C-i(ri,/5)-/C"^(ri,0)]M-i[/C-^(/3,r2)-/C~^(0,r2)] (56) 

where /C~^(r, r') denotes the kernel of the inverse of /C. It is clear that is 
invariant under a translation in both time arguments. So is U'^^\ but one has 
to carry on the calculation in (j56l) to see that explicitly; indeed, assuming 
that both time arguments lie between and /3, we have: 

rrfflW N 1 / N i(e-^^ + e'^^) z(e-^^-e'^^) \ 

where r = Ti — T2. 

Thus, the thermal propagator only depends on r, the difference between 
the time arguments, and it is naturally defined on [— 

One is usually interested in the (qq) correlation function, which here we 
can immediately read off the general expressions above, since it corresponds 
to the 11 matrix element. For < r < /5, we may derive the simpler 
expression: 

G'ii(r) = ^[(1 + nB{w))e-^^ + nB{w)e'^^] , (58) 

which is the correct result 0. This is the unique solution to the differential 
equation: 

(-dl + cu')GniT) = S{t) (59) 

subject to the condition Gii(r — P) = G'ii(r). Moreover, we can understand 
that, in our construction, the zero-temperature part appears naturally as 
the propagator for the unconstrained system. On the other hand, U^^^ is a 
solution to the homogeneous version of the equation above, and emerges as 
due to the boundary conditions. 

So far, we have seen that the previous form correctly reproduces the free 
2-point function, for all the phase space variables, when both time arguments 
are inside the [0, (3] interval. A somewhat lengthy calculation shows that G 
vanishes when one of its arguments is outside of that interval while the other 
is inside, and that it coincides with the T = function, G^'^\ when both 
are outside. Besides, in the last case, the [0, P] interval is bypassed. More 
explicitly: 



G{n,T2) = 








if 


Tl>P 


and 


0<T2<P 


GW(ri, 




if 


ri>P 


and 


T2>(3 


G(o)(ri, 




if 


ri<P 


and 


r2<P 


GW(ri- 


P,r2) 


if 


ri>P 


and 


T2<0 



(60) 



^See, for example, expression (2.32) on page 23 of [T2] . 
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and the remaining cases follow from Bose symmetry. 

This is an important property, since it allows us to perform a consistency 
check: indeed, coming back to the interacting case we considered in the 
previous subsection, one could have decided to perform the integrals in the 
opposite order, namely, first over the auxiliary fields and afterwards over 
Qa- Then, the perturbative expansion would have involved the evaluation of 
Gaussian averages of the integration term, via Wick's theorem, with G as the 
fundamental contraction. Since the interaction term involves an infinite time 
interval, one produces not only the perturbative corrections to the partition 
function (when both times are inside [0, but also a contribution that takes 
care of the normal order of the Hamiltonian (namely, when both arguments 
are outside of [0, [3\. 

3 Scalar field 

The extension of the harmonic oscillator results to the QFT of a real scalar 
field in 6?+ 1 (Euclidean) dimensions is quite straightforward. Let ip{x) = 
(^(r, x) where x = (r, x) G M*^'^"'"^), r G M and x G M.^'^\ Proceeding along the 
same lines as for the case of a single degree of freedom, we consider first the 
free case. 



3.1 Free partition function 



The free Euclidean action in terms of the phase-space variables Sq, is in this 
case given by: 



d+1 



with 



X 

1 

2 L 



tt" + |V(/^|2 



2 2 

m ip 



(61) 



(62) 



We then have to implement the periodic boundary conditions both for 
(P{t,x.) and its canonical momentum vr(r, x) 



V9(/3,x) = V5(0,x) , 7r(/?,x) = 7r(0,x) , Vx G 



(63) 



which requires the introduction of two time-independent Lagrange multiplier 
fields: ^a(x), a = 1, 2. Defining a two-component field $ = ($a), a = 1, 2, 
such that $! = (/? and $2 = ^i", an analogous procedure to the one followed 
for the harmonic oscillator yields, for the free partition function 2q(/3): 



Zo{/3) = M 



-1 



(64) 
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where ja{x) = ^a(x) [5{t - 13) - 5(r)] and: 



^ = ( f ) ' (65) 



where we have introduced h = V— + m^, the first-quantized energy oper- 
ator for massive scalar particles. Performing the integral over $, yields the 
partition function in terms of the Lagrange multipliers: 

ZoiP) = y"i)^e-5/'^'-/'^'j'54x)(x|A4.|y)Uy) ^ (66) 
with M = n{0+) + fi(0„) - - and 

n{r) ^ ( ] e-'^l-l . (67) 

Then, 



where 



-|sgn(r) 

M = ( ° ) (ns + 1)-^ , (68) 

ub = . 69 

Coming back to the expression for Zq{P), we see that: 

+ 6(x) (x|/i(nB + i)-V)6(y)]} . (70) 

By a simple field redefinition, we see that: 

Zo(/3) = det [fiB + 1) (71) 
which can be evaluated in the basis of eigenstates of momentum to yield: 

MP) = n K(^k) + 1] (72) 



where E'k = y^^+rn?. The free-energy density, Fq{[3), is of course: 
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In the classical, high-temperature limit, the path integral for the partition 
function becomes: 

Zo{P) ^ jv^e-^^^^\ (74) 

where: 

if(0 = \j d^xie^^) + |V6(x)P + m2e2(x)] . (75) 

This is, again, the usual classical expression for the partition function, with 
the Lagrange multipliers playing the role of phase space variables, and the 
integration measure being the corresponding Liouville measure. Besides, it is 
clear that the representation ( iTOl) always involves static fields, unlike in the 
Matsubara formalism. The price to pay for this 'dimensional reduction' is 
that the resulting 'action' (the exponent of the functional to be integrated) 
is spatially non local. It becomes local only in the high-temperature limit. 

3.2 Quadratic approximation 

Let us study here a simple yet illuminating example where this approach 
allows one to use non-perturbative T = information about a system as 
input for the finite temperature partition function in a quite simple way. So 
we assume that we know the exact 2-point function yV^b ^^'^o temperature. 
To be precise, since the Hamiltonian is quadratic in the canonical momentum, 

('2) 

we only need to know , since it is possible to show that when one of the 
fields is replaced by the canonical momentum, the result is multiplied by the 
corresponding frequency (in Fourier space). Indeed, this can be shown, for 
example, by performing the exact integral over the canonical momentum. 
Then the effective Hamiltonian corresponding to this term has the form: 

HefAi) = \ I [6(x)Ci(x - y)ei(y) + 6(x)C2(x - y)6(y)l (76) 

where the Fourier transforms of the coefficients Ci and C2 are 

CM = |^(l-e-^^^-o)VVn(fco,k) 

C'2(k) = ^ / ^ (1 - e-'^'^") kl m,{ko, k) . (77) 

Note that these two coefficients that determine the contribution of the quadratic 
term to the partition function, could be rather involved functions of the func- 
tions of a coupling constant, since we are not perturbing but just assuming 
that we only consider the information contained in the full propagator. 
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Assuming now that the exact propagator is such that only one particle, 
with energy E{k), exists, one knows that there can only be poles associated to 
them in Wn. Moreover, we assume that, as usual, renormalization conditions 
have been imposed such that the residue is 1. Then we obtain: 

C.(k) = (1 _ , (78) 

and the partition function becomes: 

^-(/^) = n i_e-mk) ^ (79) 

which is of course the one of an ideal Bose gas. But the important difference is 
that one is obtaining it by putting the information contained in the knowledge 
of the non-perturbative spectrum. It could be the case, for example, of a 
model where the mass is generated by a non-perturbative mechanism. Or one 
could have a theory with more than one pole, corresponding for example to 
different bound states. These non-trivial poles, obtained in the T = theory 
are then directly taken into account by this contribution, as the first terms in 
an expansion in powers of the auxiliary fields. Of course, one has to include, 
in general, also non-quadratic terms. The quadratic approximation could be 
justified, for example, within the context of a large-A^ approximation. 



3.3 Perturbation theory 

We study here the alternative representation for the partition function in an 
interacting theory. For the sake of clarity, we present this topic within the 
context of a concrete example: the real scalar field in 3 + 1 dimensions, with 
a self-interaction of the quartic type. The action is then 

S = So + Sj , 5, = ^ y d^x^\x) . (80) 

As in the standard formulation, we want to calculate Z{i3) in a power 
series of the coupling constant, A. It should be clear that the proper way 
to do that here is to calculate W[J] to the desired order, and from there 
to obtain the corresponding effective Hamiltonian. The resulting expression 
for the partition function as a functional integral over the auxiliary fields 
can then also be expanded (up to the same order W[J] was calculated). In 
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the course of such a program one has to face the issue of the UV infinites 
that pop up in the calculation of loop diagrams for the (zero-temperature) 
object W[J]. We shall assume that all those infinities are renormalized in 
the usual way at zero temperature. In particular, our T = vertices are 
normal ordered, so that the zero-temperature tadpoles are, to begin with, 
absent. The finite-temperature tadpoles will, nevertheless, appear in the 
alternative description, as we shall see. Namely, the functional integral in 
terms of the auxiliary fields is UV finite. In this way, the thermal corrections 
are completely disentangled from the issue of renormalization. 

We now study the problem of calculating the effective Hamiltonian in 
terms of the perturbative expansion for W[J]. As in the case of one degree 
of freedom, the effective Hamiltonian is related to W[ J] by: 

HeffiO = W[zj(x)] (81) 

where VV[J] is the generating functional of connected correlation functions 
$a, at T = 0. Moreover, we shall assume that, in the previous expression, 
>V[0] = 0, since any vacuum contributions would be cancelled by the nor- 
malization factor J\f. 

The perturbative expansion of W will be denoted by W = >V^°^ + yV^^\ 
with 

Wi[J] = W^^V] + W(2)[J] + ... (82) 

where the index denotes the order in A of the corresponding term. This 
yields the corresponding expansion for the effective Hamiltonian, H^jf = 
uf^j + H^ff-, and one can then find corrections to the partition function, 
or the free energy F = — ^InZ, by evaluating the corresponding Gaussian 
averages. Indeed, 

F{P) = fW(/3) + F(^)(/3) (83) 

where 

F(^)(/3) = -lln(e-^</(«)) (84) 
where the average symbol is defined by the quadratic weight: 

To fix ideas, we do that first for the simplest non-trivial order, i.e., A. 
Then one has the first-order contribution to F: 

F^'HP) = (<^/(0) , (86) 
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where the Gaussian average requires the knowledge of the elementary aver- 
ages which involve two auxiliary fields. The only non-trivial ones are: 



(6(x)ei(y)) 
(6(x)6(y)) 



^ik-(x-y)_ 



,ik-(x-y; 



1 _ e-/3<^(k) 
1 



-/3^^(k)^ 



(87) 



(27r)3 cj(k)(l-e- 

On the other hand, a standard T = calculation shows that the first- 
order term in the expansion of W is: 

X r . r r . ^ i4 



4! 



d X 



d'yY.^f!{x-y)Uy) 



where (x — y) is the free T = propagator for the real scalar field and 
its canonical momentum, namely 

G^'^{x-y) = (89) 

which may be conveniently represented in terms of its Fourier transforms 
G^^^{k), in a matrix representation where (1 corresponds to ip and 2 to tt): 



G^'^\k) 



J^2 _|_ ^2 



1 



-iko kf 



iko 

2 




(90) 



Using now the rule that maps terms in the expansion for W into like ones 
for Heff, we see that: 



= ^ jd^x[jdSY.GTAx-y) 



A 



.(0), 



X ea(y)(5(yo-/9)-5(yo)) 



(91) 



Then, using tildes to denote the Fourier transforms of the auxiliary fields, we 
see that each one of the factors that appear integrated over x above, can be 
put in the following form: 

j d'y J2 Gf] {x - y)Uy) Hvo -13)- 5(yo)] 



d^k 1 



g*k-x r -a;(k)|xo-/3| _ ^-w(k)|a;o| 



(27r)3 2uj{k) 
d^k 



]ei(k) 



^ik x " 



(2vr) 



[cr(xo)( 



-a;{k)|xo| 



a(a;o-/?)e--('')l^«-^l]6(k) , (92) 
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where a{x) denotes the sign of x. 

Introducing the expression above into (lOTl) we find, in a natural extension 
of the notation used for the case of one degree of freedom: 

= 5/[n||r](2')'*'''(Ek.)wiS.(k,,...,k,) 

' ^ i=l ^ ' i=l 

X 4,(ki)...ea4(k4) (93) 

where HS^'^^ is the kernel for a quartic term; the '1' has been written to 
pinpoint the fact that it has been calculated to the first order. The explicit 
form of the kernel elements depends of course on the indices considered. 
Modulo permutations, the only inequivalent possibilities are summarized in 
the following results: 



^ 4 , 4 4 -il3ui 1 

<';(k.k.k„M ^ A / (n^) 2.*(E^.) n [^Ti^^jj 

1=1 7 = 1 1=1 ' ^ ' 



j = l 1=1 



'^iii2(ki,k2,k3,k4) 



i=l 



X 



n 



^^[6-^^"^ - 1) 



(94) 



(95) 



4 4 

7^itS(k„k„k3,k4) = A /(n^)27r5(x:^,) 

i=l 7 = 1 



X 



n 



n 



VIV. 



(96) 



M222(ki,k2,k3,k4) 



i=l 



n 



pz/,.(e-'^'^'- - 1) 



vl + uj\\,,) AIL z/2 + ^2(k 



and 



(97) 



«s,k..k.k3,k., ^ A / (n|^)2..E.,, n[^^fc.^ 

j = l i=l « ' V 



i=l 



(98) 
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Upon application of Wick's theorem for the calculation of the average of the 
first order effective Hamiltonian, we see that only the terms with an even 
number of legs for each field yields a non- vanishing contraction; namely, only 
the terms Tin'n? ^1122 ? ^2222 enter into the calculation. Besides, the first 
and third of these carry a factor of 3 because of the number of unequiva- 
lent contractions, while for the second one there is a 6 due to the different 
permutations of the (different) indices. 

Using the explicit form of the contractions, and integrating over the fre- 
quencies, one sees, after a somewhat lengthy, but nevertheless straightforward 
calculation, that the proper result is obtained. Namely, 

where V is the spatial volume of the system. 



3.4 Generating functional 



Including a source in equation ( 1641 ) it is straightforward to obtain the gener- 
ating functional for the scalar field in d + 1 dimensions in a analogous form 
than for + 1 dimensional theory. Working in Fourier space for the spa- 
tial coordinates, we see that the entire + 1 dimensional procedure applies. 
Then Gab{x,y), the correlation function for the scalar field and its canonical 
momentum becomes: 

G^,{x,y) = gS)(x -y)- ui^{x - y) , (100) 

where G^^ij{x — y) is the zero temperature correlation function, whereas 
{x — ?/) is a temperature dependent function. 
The explicit form of these functions is: 



G^h- -y) = (101) 



where: 



and: 



.(0). , _ e--'^^ / 1 lujk 



which are exactly equal to the + 1 dimensional case equations on which we 
have to replace w hy cuy^ = v^FT^r?. 
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4 Dirac field 



The final example we consider is a massive Dirac field in c? + 1 spacetime 
dimensions. The procedure will be essentially the same as for the real scalar 
field, once the relevant kinematical diff'erences are taken into account. The ac- 
tion 5*0 for the free case is given by = / d'^~^^xip{^ + m)'ilj where ^ = 'j^dfj^, 
ll = li^ and {7^, 7^} = 26^^. 

We then impose antiperiodic conditions for both fields: 

^(Ax) = -z^(0,x) , ^(/3,x) = -^(0,x) (104) 

as constraints on the Grassmann fields. Those conditions lead to the intro- 
duction of the two 5— functions: 

Z^iP) = y I)V^P^5(^(/5,x)+V^(0,x)) 5(^(Ax)+^(0,x)) 

X exp -S^iij,^}) . (105) 

Since the Dirac action is of the first-order, the introduction of two constraints, 
and two Lagrange multipliers, appears in an even more natural way than for 
the previous case. Those auxiliary fields, denoted by x(x) and x(x) must be 
time-independent Grassmann spinors. The resulting expression for is 
then 



where r] and r] are (Grassmann) sources depending on x and x through the 
relations: 

V{x) = x(x) [6ir - /?) + Sir)] , r/(x) = x(x) [5(r - P) + (5(r)] . (107) 
Integrating out ijj^ip, we arrive to: 

Z/,iP) = J VxVxexp[-PH,ff{x,x)] (108) 



where 



with: 



Heff{x,x) = jd'^ j d''yx{^)H^'\^,y)x{y) (109) 



if(2)(x,y) = (x,0|(^ + m)"V,0) + (x,/3|(^ + m)-i|y,/3) 
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+ (x, 0|(^ + m)-i|y, P) + (x, + m)-'\y, 0) 



P 



25/(0,x-y) +5/(/?,x-y) +5/(-/3,x-y) 



(110) 



On the last line, Sf, denotes the Dirac propagator. A quite straightforward 
calculation shows that 



H{^,y)=-{^\uil-nF)-'\y) 



(111) 



1 + 6' 



/3n 



is the Fermi-Dirac distribution function, written 



where fip = 

in terms of h, the energy operator (defined identically to its real scalar field 
counterpart); m is a unitary operator, defined as 



h 



u = 



D 



hf) = J ■ V + m . 



Then we verify that: 

Zo^(/3) = detM det-^[(l -n^)l] , 
(I = identity matrix in the representation of Dirac's algebra) 



(112) 



(113) 



(114) 



with E{p) = y/p^ + and = dimension of the representation (we have 
used the fact that detu = 1). 

Again, the procedure has produced the right result for the partition func- 
tion, with a normal-ordered Hamiltonian. On the other hand, for a Dirac 
field in a static external background corresponding to a minimally coupled 
Abelian gauge field the Aq = gauge, we have 



X 



ip{x) ( ^ + 267- A(x) + rnjip{x) 



(115) 



The assumed r— independence allows us to carry on the derivation described 
for the free case, with minor changes, arriving to the expression: 



Z^{(3) = detM(A)det-^(ni.(A)l) 
= e*^(^)det [(l + e-^'^(^))l 



(116) 



23 



where 
and: 



h{A) = V-D2 + m2 , T> = V-ieA 
det (7 ■ D + m) 



JK(A) 



dety^D2" 



(117) 



(118) 



Notice that the factor det 



can be formally diagonalized in 



(l + e~MA))i 

terms of the energies Ex{A) in the presence of the external field. Thus we 
arrive to the expression: 




-pEx{A) 



(119) 



The factor e''^^-^\ on the other hand, is topological in origin, as it depends 
on the phase, K{A), of the determinant of hf). On the other hand, fiD 
may be regarded as a kinetic operator in one fewer dimension. For Dirac 
fermions, we know that the phase of det ho can be non-trivial only when d 
is odd, i.e., when d+1 is even. However, the 7-matrices appearing in det ho 
form a reducible representation of the Dirac algebra in d dimensions, with 
the matrix 7,- relating every eigenvalue to its complex conjugate. Thus, as a 
result, the phase K{A) vanishes. Of course, a non-vanishing result may be 
obtained for other fermionic systems, like Weyl fermions for d + 1 = even. 



5 Conclusions 

We have shown that, by introducing the periodicity conditions as constraints 
for the paths in the Euclidean path integral for the T = vacuum functional, 
one can obtain a novel representation for the partition function. These con- 
straints should be applied on fields and canonical momenta, and when they 
are represented by means of auxiliary fields, they lead to an alternative, 'dual' 
representation for the corresponding thermal observable. 

Since both phase space variables should be constrained, one has to work 
in a first-order formulation; this is automatically satisfied in the case of a 
Dirac field, but it requires a little bit of care in the case of the real scalar 
field. 

The resulting representation for the partition function may be thought of 
as a dimensionally reduced path integral over phase space, similar to the one 
of a classical thermal field theory, with the auxiliary fields playing the role of 
canonical variables, but with an effective Hamiltonian, -ffe//, which reduces 
to the classical one in the corresponding (high-temperature) limit. 
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We analyzed the main properties of this representation for the cases of 
the real scalar and Dirac fields, two typical examples that have been chosen 
for the sake of simplicity. It is not difficult to generalize the representation 
to the case of systems containing fermions interacting with bosons. For 
example, assuming that S{ijj, ip, $) is the first order action corresponding to 
a real scalar interacting with a Dirac field, we define the T = generating 
functional W(C, C, J) by: 

W{CX,J) = ^ j V^V^jV^ 6-^^^^'^'''^+' '''^"''i^'^^^^^-^-'^-') , (120) 

where the sources are arbitrary. Then, -ffef/ can be obtained from the ex- 
pression: 

^^e//(x,X,0 = W(r/,r/,^j) , (121) 

where ?7, fj and J are (the already defined) functions of the Lagrange multi- 
plier fields X, X ctnd ^. 

We have shown how the effective Hamiltonian can be constructed by as- 
suming the knowledge of the corresponding T = generating functional of 
connected correlation functions. If this knowledge is perturbative, one recov- 
ers the perturbative expansion for the thermal partition function. However, 
the most important applications of this formalism are to be found in the case 
of having non-perturbative information about the T = correlation func- 
tions: here, it is quite straightforward to incorporate that knowledge into 
the formalism, and to compute thermal corrections from it. 
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